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Abstract 

Effective Liouville operators of the first- and the second-order symplectic in- 
tegrators are obtained for the one-dimensional harmonic-oscillator system. The 
operators are defined only when the time step is less than two. Absolute values 
of the coordinate and the momentum monotonically increase for large time steps. 
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Symplectic integrators are widely used not only for the Hamiltonian dynamics. 
Higher-order ones are also studied as the product of exponential operators [1, 2] and 
are performed in various areas. However, the upper limit of the time step where the 
scheme works is not studied yet. Suzuki [3] studied the convergence of decompositions 
of exponential operators in a Banach space. Bourbaki treated only for the complete and 
normed Lie algebra. In the present Letter, we discuss the convergence of the scheme for 
a specific unbounded system. 

Let us consider the following first-order symplectic integrator for the one-dimensional 
harmonic oscillator described by the Hamiltonian Ti = {p^ + 9^)/2: 

f p' =p-xq, 
\ q' = Q + xp' ■ 

No anomalies seem to exist because the above discretization scheme is defined for any 
time step x. However, the absolute values of p and q monotonically increase for a; > 2. 

The above anomaly is explained by the divergence of the "effective" Liouville oper- 
ator which corresponds to the symplectic integrator. The discretization scheme Eq. (1) 
is expressed by 

= (1 + xA)(\ + xfi) f M = exp(xA) exp(x5) { ^ j 
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with two matrices 



10' 



because of = 0, = 0. The equations ABA = -A and BAB = -B hold, and then 
we have 

exp{xA) exp{xB) = cxp (^A + B + ^[A, B]^ | , (2) 

where the function F{x) is defined by 



(n!) 



n2 



X 



2n 



with the convergence radius two [4]. 
Defining an operator Li by 

L,^A + B + ^[A,B]=(^\ 

we have iF{x)Li as an effective Liouville operator which corresponds to the symplectic 
integrator Eq. (1). When a conserved quantity Ei is expressed as a quadratic form of p 
and q using a real symmetry matrix Mi such as 

E, = {p q)M, ^ j , 

the matrix MiLi is alternate, because the time derivative of Ei is obtained from 

^ = (p g)%F(a;)Mi ^ ) + (p g)M,F(a;)Li ^ ) 



Therefore, with 
we have 

Ml 

that is. 



MiLi = idet(Li) ( 5 ^ 







1 j 1, -1 f ; 2 ^ -f 1 



El = ^{p^ - xpq + q^). 
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Within the convergence radius, the above Liouville operator iF(x)Li describes con- 
tinuous dynamics along an clUpsc in the (p, g)-plane with a constant value of Ei. Di- 
vergence of F{x) can be regarded as a break down of the continuous movement. 

Here we derive Eq. (2). 

Expanding log(exp(a;y4) exp{xB)) as the formal power series of A and B, the following 
four kinds of terms appear: 

ABA ■■■BA, 

2n+l 

BAB---AB, 

2n+l 

ABA---AB, 

2n+2 

BAB ■■■BA 

^ V ' 

2n+2 

because of = 0, B"^ — 0. Coefficients of the above terms are given by Goldberg's 
theorem [5] as follows: 

(-1)" 



(2n + 1)! 



2n+l 



' (2n+l)! 

n!(n + l)! _ ( ^^n ^2n+2 

> (2n + 2)! ^ - y 2(2n + l)! 

\2 



^ ' (2n + 2)! ^ ' 2(2n + l)! 

Because of ABA = -A, BAB = -B, we obtain 

ABA---BA = {-If A, 

2n+l 

BAB ■ ■ ■ AB^ = 

2n+l 

ABA^^_AB = (-l)MS, 

2n+2 

• • • = {-lyBA, 

2n+2 

and we arrive at 

log(exp(x>l) exp{xB)) 
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^'(2n + l)! (2n + l)! 2(2n+l)! 2(2n + l)! J 



n=0 

oo 



= xF(x)(^A + B+^lA,B]y 

Note that log{exp{xA) exp{xB)) is expressed only with A, B, and [^4, B] because of 

[A,[A,B]]^2A, [B,[A,B]]^-2B. 

Next we consider the second-order symplectic integrator, which is expressed with the 
two matrices A and B as 

= exp exp(x>l) exp ( ^ ) . (3) 

Let us define by 

S{x) = exp exp(x^) exp , 

we obtain S{—x) = S{x)~^ because of S{x)S{—x) = 1. When we express S{x) as 

S{x) = exp/(x), 

the following equation holds: 

exp/(-a;) = S{-x) = S{x)~^ = exp{-/(x)}, 

which indicates that f{x) is an odd function. Therefore, when we expand f{x) as the 
formal power series of A and B, the following two kinds of terms appear: 

ABA---BA , (4) 

2n+l 

BAByAB^ . (5) 

2n+l 

Because of ABA — —A, BAB = —B, we can express f{x) as 

f{x) = x{F^{x)A + F2{x)B) 
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with two even functions Fi{x) and F2{x). 

When we expand log(exp Xi exp X2 exp X3) as the formal power series of Xi, X2, and 
X3, coefficients of the terms are given by an extended version of Goldberg's theorem 
[6, 7]. Because we later replace Xi, X2, and X3 as 

Xi^ X^ = -^B, 
X2 = xA, 

two kinds of terms 

-'^2-'^n-'^2-'^i2 ^2-^^13 -^^2 • ■ ■ -'^2-'^j„-'^2, (6) 

^ V ' 

2n+l 

^11^2X12^2X13X2 ■ • • -'^2-'^i„-'^2-'^i„+i (7) 
^ V ' 

2n+l 

contribute to the terms (4) and (5), respectively, with ij = 1 or 3. According to 
Goldberg's theorem, coefficients do not depend whether Xi^ between two X2^s is Xi 
or X3. Then the coefficient of the term (6) is given by 



(2n + l)! 

for all Xy. The coefficient of the term (7) is also given by 



(2n+ 1)! 



for = (1,1), (3,3), 



{n-l)\in+l)\ 



(2n+l)! ^""^ (ii,in+i) = (l,3), (3,1), 

if n > 1. Replacing Xi, X2, and X3 with A and B, and using ABA — —A and 
SAB — —B, we arrive at 

that is, 

exp (^b\ exp(a;A) exp (^b\ = exp \xF{x) (a+(i- ^] B] \ . 



5 



As for the first-order scheme, defining L2 by 

we can regard iF{x)L2 as an effective Liouville operator which corresponds to the second- 
order symplectic integrator Eq. (3). Then we obtain the conserved quantity E2 as 

Finally we note that the exponential function absolutely converges and the product of 
exponentials always exists, which corresponds to the fact that the discretization scheme 
as Eqs. (1) and (3) can be defined for all time steps. However, the absolute values of p 
and q monotonically increase for large time steps although Jacobian of the time-evolution 
operator is unity 
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